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In these lectures we look for parallels between symmetry breaking in the early universe and 
condensed matter systems, and discuss experiments that display these. 



I. OVERVIEW 

Our understanding of the early universe is strongly 
hampered by the difficulty in making clean inferences 
about the nature of its underlying field theory from the 
observational data. The goal of the COSLAB programme 
is to find general physical principles common to particle 
astrophysics (Cosmology) and condensed matter physics 
(the Laboratory) which permit direct experimental test- 
ing in the latter that is impossible in the former. 

Phase transitions are ideal in this regard, occurring 
in both condensed matter systems and the early universe 
and, in the former case, being directly testable in the lab- 
oratory. In these talks I shall examine some of the paral- 
lels that phase changes suggest, their consequences and 
their experimental confirmation. Because of the breadth 
of these talks, in many cases I have cited review articles, 
rather than the original references. Brevity requires oc- 
casional over-simplification. I apologise to authors who 
feel unfairly omitted. In particular, I refer the reader to 
the proceedings of earlier COSLAB meetings PIQ- 



II. PHASE TRANSITIONS 

Phase transitions are changes of state, associated with 
changes of (internal) symmetries of the system. Although 
they occur in both the early universe and condensed mat- 
ter systems, their status is different in the two cases: 

In the case of the early universe, we have many transi- 
tions anticipated on general grounds. However, although 
we have little knowledge of the symmetries at very high 
energy scales, we can make some generic predictions. For 
example, simple models typically produce monopoles and 
cosmic strings, both of which are observable. Unfortu- 
nately our understanding of the non-equilibrium dynam- 
ics of the very early universe that would enable us to 
predict the number of each is very poor (but see Paul 
Shcllard's lectures in these proceedings). 

On the other hand, for condensed matter, many transi- 
tions are known, for which our detailed knowledge is very 
good. However, there has been limited experimentation 
on them in directions parallel to the interests of the cos- 
mologists, which emphasise their non-equilibrium nature. 
For these reasons, phase transitions in condensed matter 
provide an ideal testing ground to see to what extent one 
community can inform the other. 

In condensed matter systems familiar changes of state 
include: 



• Non-Magnetic Magnetic 

• Normal Fluids Superfmid 

• Normal Conductance Superconductance 

• Bose-Einstein Condensation 

• Liquid crystal formation 

There is no hierarchy within these transitions. In 
particular, they are true transitions, rather than rapid, 
but smooth, cross-overs. The situation is very different 
for the early universe, for which early times imply higher 
energy which, in turn, implies greater symmetry. In 
consequence, a cooling universe generates a hierarchy of 
transitions: 
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In practice, greater symmetry implies less knowl- 
edge (the more stars the better). The 4-star hadronic 
transition is best known. Its cold phase (the hadrons of 
accelerator physics) is extremely well understood. The 
transition is unique among early universe transitions in 
that the hot phase, the quark-gluon plasma (QGP), can 
be simulated experimentally, in some sense, by colliding 
heavy ions at high enough energy. Our theoretical 
understanding of deconfinement is also good, even if we 
cannot pin down all the details. However, there are also 
high baryon number superconductor/superfluid transi- 
tions that could be important in neutron stars, which 
we cannot recreate. As we shall see, the transitions are 
often cross-overs, but this will turn out to be one of the 
less important differences. 

For the 3-star electroweak transition particle accelera- 
tors have given us the main properties of the cold phase 
(but for the Higgs sector). The hot phase has been esti- 
mated theoretically to suggest a complicated transition, 
again a crossover, but in this case there is no direct ex- 
perimental confirmation. 

For our 2-star GUT transitions there is no direct ex- 
perimental evidence on either side of the transition tern- 



2 



perature(s). Our belief in their likelihood is largely moti- 
vated by the observed convergence of coupling strengths 
at GUT scales due to the screening behaviour of field fluc- 
tuations. Our 1-star supersymmetry transition, which 
balances fermionic against bosonic degrees of freedom, 
has even less direct support. We are motivated by the so- 
lution that supersymmetry brings to the hierarchy prob- 
lem (the need to keep GUT and electroweak symmmetry 
breaking scales separate) and the fact that theories that 
incorporate gravity at a quantum level require supersym- 
metry. At the Planck scale we know nothing reliably, but 
it sets the absolute scale against which all others are com- 
pared. 



III. PATTERNS OF SYMMETRY BREAKING 

For true transitions a change of phase is typically iden- 
tified by the vanishing of an order parameter (or of one 
or more components of an order parameter vector or ma- 
trix). We have two cases to consider: 

• Bosonic Theories: The order parameter fields 
are bosonic 4> a {x) from which we construct or- 
der parameters as ground state expectation val- 
ues 4> a = (4>a)- Typically the change in symmetry 
is implemented by symmetry breaking in a scalar 
(semi)classical field potential. 

• Fermionic Theories: Fermionic fields have zero 
ground state expectation values. The order param- 
eters are constructed from expectations of bosonic 
fermion bilinears, such as (ipip) (quark condensate) 
or (Cooper Pair). Often we can write the 
underlying fermionic theory as an effective bosonic 
order-parameter field theory, so that we can still 
use simple scalar potentials to characterise phases. 

There are no problems for global symmetry transforma- 
tions, exemplified by 4>(x) — > e lea <j){x), where a is con- 
stant in space-time. However, for local symmetries, most 
simply of the form <p(x) — > e lea ( x ' 4>{x) , with space-time 
dependent a(x), the symmetry requires the existence of 
gauge fields that undergo their own transformations 

For global symmetries the exemplary groups are the 
Orthogonal Groups O(N), with a vector 4>{x) of order 
parameter fields <p a (x), a = 1, 2, ...N, transforming under 
0(N) as <j)(x) — ► R<j)(x). R is a real N x N orthogonal 
matrix, with RR = I. If detR — 1 the symmetry group 
is SO(N), rather than O(N). The simplest invariant 
from which to construct the potential is (f> 2 . [We should 
also consider the additive Groups Zn, for which (p(x) — > 
U(j)(x); U N = I. The most common of these is Zi : 
(f)(x) — > —(j>(x) for a single real field.] 

For local symmetries the exemplary groups are the 
Unitary Groups U (N) , although some global symmetries 
are of this form (e.g. our U(l) example above). The or- 
der parameter fields most simply transform under the N- 
dimensional vector representation as (j>(x) — » U(x)<fi(x), 



where WU — I. With the iV-component <f> complex, 
the basic invariant is (jy(x)(p(x). If detU = 1 then the 
group is SU(N). As we noted, the requirement of lo- 
cal symmetry invariance demands the existence of gauge 
fields A a ^{x), which transform under the adjoint repre- 
sentation of the group. If we restrict ourselves to uni- 
tary groups U = U{x) there are N 2 gauge fields A% (x) 
(a, 6=1,2, N) before gauge fixing. For SU(N), with 
A% (x) traceless, there are N 2 — 1 gauge fields. 

Note that [7(1) = SO (2) : $<j> = cj>l + <f%. Also note 
that SU(2) « 50(3), 50(3) x 50(3) sa 50(4), where 
the « means identical behaviour for infinitesimal trans- 
formations, but different large-scale behaviour. Other 
groups (Sp(N), G2, Ei, ...) are rarely used. 



A. General Symmetry Breaking 

Transitions occur when the system has degenerate 
ground states. That is, the symmetry group G of the 
Hamiltonian, or the Action (Lagrangian), is not the sym- 
metry group H of the ground states of the system. We 
write G H . H is determined by enumerating the 
ground states cf>Q. For all h S H C G, we have h(po = 4>a- 
This enables us to identify the manifold of ground states 
M as the left cosets M — G/H. 

From our comments above, we shall restrict ourselves 
to spontaneous symmetry breaking (SSB) visible in a 
classical, or effective, bosonic potential. However, some- 
times the symmetry of the Hamiltonian (and the degener- 
acy) is only approximate, and we may need additional ex- 
plicit symmetry breaking by mass terms, etc. Such terms 
can induce significant quantitative changes. In particu- 
lar, true transitions are replaced by crossovers. Details 
on symmetry breaking can be found in many sources, 
particularly Kibble in [J, Q and we can only recreate the 
rudiments here. 



B. Some simple examples 

The simplest symmetry breaking is that of the Zi (</> — > 
—(f>) symmetry of the generalised Ising model of a single 
real field 0, with double-well potential V{<f>) — \\{<t> 2 — 
ij 2 ) 2 , and action 



5 = I d 4 x 
or free energy 

F= I d 
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The maximum of V at <p = is unstable, but on expand- 
ing about the minimum <f> = rj, as <p = r\ + h, we see 
that V(4>) — \rr? h \{ 2 + (non- quadratic terms) describes 
/vfield excitations of mass mh — r]V2X. Expanding about 
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= — 7/ gives an identical spectrum. Z2 symmetries have 
a natural role in Grand Unified Theories (GUTs) of the 
early universe, where they provide a simple selection rule 
against the rapid decay of the proton. 

The next model, to which we shall return repeatedly, is 
the Goldstone model of a complex field = 01+702, with 
Mexican-hat/wine-bottle potential V(<f>) = |(|0| 2 - r) 2 ) 2 
(and kinetic/gradient terms changed correspondingly in 
Q and J5J). The global 0(2) symmetry of rotations in 
the complex plane is totally broken by the inequivalent, 
but equally acceptable, degenerate ground states with 
|0| = 77, which form the manifold M. = S 1 (the 1-sphere, 
or circle). Expanding about — 77 say, shows that the 
particle spectrum is that of a massive (Higgs) particle, 
mass rnt as above, corresponding to radial oscillations in 
the complex 0-plane, and a massless (gapless) mode, the 
Goldstone particle, corresponding to translations of the 
field along S 1 . 



The natural generalisation of the Goldstone model is 
to N scalar fields (j> a , a = 1,2,...,N, with potential 
V = |(^-r? 2 ) 2 = |(0 2 - 77 2 ) 2 , where # = # + $ + 
The O(N) symmetry is broken as O(N) -> 0(N - 1), 

with a manifold of groundstates 2 = r/ 2 , the N — 1 
dimensional sphere M — S ,JV_1 . If we expand about 
01 = r), 02 = 03 = ••• = 0, the oscillations describe 1 mas- 
sive (Higgs) field in the radial field direction and N — 1 
massless (Goldstone) fields along the surface of the sphere 
of groundstates. Such a model is termed the O(N) Lin- 
ear Sigma Model (LaM). The O(N) Non-Linear Sigma 
Model (NLaM) is an extreme version of this, describ- 
ing classical rotors, in which, by fixing 2 — r/ 2 , we 
eliminate the Higgs field as a degree of freedom. There 
are no known Goldstone modes in particle physics, al- 
though there are very light degrees of freedom (e.g. pi- 
ons), but Goldstone modes arise naturally in condensed 
matter physics as phonons. 

For gauge theories it is not sufficient to just examine 
the potential V(0), because the local nature of the gauge 
transformations enforces additional interactions through 
the covariant derivatives that extend the kinetic terms. 
The simplest case is the U(l) — > 1 Abclian Higgs model 
for complex field and gauge field A^. The action is 



S = 



d D , 



\D„ 



I 2 - v 2 ) 2 7 V 4 " 



where = g§- 



ieA^ is the covariant derivative, and 



F^ u the electromagnetic field tensor. Although M. = S 1 , 
as for the Goldstone model, there are no Goldstone 
modes. The rearrangement of the field degrees of freedom 
gives rise to one real Higgs particle, and one massive vec- 
tor field. In some sense the Goldstone mode is swapped 
for the transverse degree of freedom of the vector field. 



IV. SYMMETRY BREAKING IN CONDENSED 
MATTER 

The simple forms of SSB given above are most eas- 
ily realised in condensed matter theory. We shall largely 
restrict ourselves to continuous transitions for which rel- 
evant experiments have been performed. 



A. Superfluid *He: 

This is the global 0(2) Goldstone model already dis- 
cussed, with complex order parameter field (f> = pe l6 . In 
the 2-fluid model n s = p 2 is the superfluid density and 
v s = ^V6* the superfluid velocity. The Goldstone mode 
describes sound. As before, M = S 1 . 



B. Low -T c Superconductors 

The bosonic effective order parameter field <f> is derived 
from the L = S = Cooper pairs of electrons (with mo- 
menta close to the Fermi surface) as <j) ~ (ipitpi)L=s=o- 
This is an Abelian Higgs model with Free Energy 
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(3) 

in the absence of external fields, where e* = 2e, 771* = 
2m e characterise the Cooper pair. The Meissner effect, 
whereby the magnetic field penetrates into the bulk su- 
perconductor, has its penetration length determined by 
the mass of the vector field in the broken phase. 



C. High- T c Superconductors 

High-T c superconductivity is a complicated phe- 
nomenon. For our purposes, we adopt an idealised ex- 
plicitly broken SO(5) model in D = 2 dimensions 0. 
The basic idea is that doping an antiferromagnet leads 
to d-wave superconductivity. The bosonic effective order 
parameters are, again, constructed from fermionic bilin- 
ears. 

We begin with a global 0(3)af LctM for antiferro- 
magnetism with order parameter field (staggered mag- 
netism) 77 with potential V(n) — j(n 2 — l) 2 . This 
is extended to a five-component order parameter field 
N = (</>i, ni, 712, 713, 02), with 0(5)-invariant potential 
V{N) = ±{N 2 - l) 2 . [We can work equivalently d H 
with a NLcrM, in which |AT| = 1.] Ultimately we shall 
couple <j> — (0i + 702) locally to the EM field as in Low- 
T c superconductors. 

In the first instance M. = S 4 . We break the 0(5) 
invariance explicitly to 0(3)^ x SO(2) by the addition 
of antiferromagnet and doping terms to the potential, 
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most simply as [j| 



V(cj),ff) = ^a s \cb\ 2 + ^a A n 2 



?2\2 



(4) 



where as,au < 0. We then couple <j> to EM as 

h 2 





" h 2 




j d 2 x 




(- iV -fc A )^ 


2m* 





2m" 



F= I d~.v - — [ -/V- —A. ! < - + — (Vriy+Vi&n) 

(5) ' 

Increasing the doping drives \ag\ > |<Xr|, making the 
U(l) ~ SO(2) superconducting direction the global min- 
imum. That is, 0(5) -> 0(3)af x 50(2) -> 50(2), with 
= 5 1 . This 17(1) - 50(2) is then broken as for the 
low-temperature superconductors . 



D. Superfluid 3 He 

3 He is a Fermi Liquid, which can become superfluid by 
the formation of p-wave (L=S=1) 'Cooper pairs' of 3 He 
atoms. L and S are uncoupled at short distances, to give 
a global symmetry group G = SO(Z)l x 50(3)s x U(1)n- 

The effective order parameters 6j form a 3 x 
3 matrix A a i(x), formed from the Fermi bilinears 
(i/j(x)iP{x))l=s=i- The label i = 1,2, 3 is the orbital 
angular momentum label and a = 1,2,3 the spin label. 
The £7(1) at describes the overall phase freedom. Above 
the transition all the elements of the matrix have zero 
values. Below the transition, some of these quantities 
become non-zero. The symmetry of the order parame- 
ter after the transition corresponds to the manifold of 
symmetries which remain unbroken. 

The free energy of these states can be expressed in the 
framework of the phenomenological Ginzburg-Landau 
theory by a potential 



V GL (A) 



~ aA *a,t A a,i 



+02A ai A at iA b jAb,j + f5$A ai A bi A a jA b j 
fcK,i A b,iKjKi + 05A* a:i A b>i A b>j A* ad (6) 



The different possible symmetries of the order parameter 
A a i are identified with local minima and saddle points in 
this 18-dimensional energy surface. There are two stable 
phases; 

• The A phase, in which 

SO(3) s xSO(3) L xU(l) N -» SO(2) Sz xU(l) L ^ N/2 xZ 2 . 

The manifold of ground states is Ma — G/Ha = 
S 2 x SO(3)/Z 2 . The order parameter in the A 
phase ground state is anisotropic in both spin and 
orbital spaces (the 'axial' state). : most simply, it 
takes the form A° a i = A a z a (ii + ijji), where x, y, z 
are unit vectors in the x, y, z directions respectively. 

• The B phase, in which the orbital and spin angular 
momenta are locked together as 



50(3) 5 x 50(3) L x U(1) N -» SO (3) 



S+L- 



The manifold of ground states is now Mb = 
G/H B = 5 1 x 50(3). As a result, A a ^ resembles a 
rotation matrix. Specifically, in the bulk B phase, 
A a ^ reduces to the arbitrary orthogonal rotation 
matrix A a>i = AR aii e l ® . 

The energy balance between the A and B phases is 
determined by the relation between the parameters /3j. 
At zero pressure, the B phase corresponds to the abso- 
lute minimum, while at pressures above 20 bar there is a 
temperature range in which the A phase is preferred. 



E. Other systems 

There are other systems whose transitions are under- 
stood well, which potentially have parallels with the early 
universe. In particular, we would cite 

• Uniaxal nematic liquid crystals, for which the or- 
der parameter in the nematic phase is the non- 
oriented director vector with ground state mani- 
fold RP 2 = S 2 /Z 2 . The transition is first order 0. 
However, when considering symmetry breaking at 
the interface of an isotropic-nematic transition the 
anchoring of the director at the interface forces it 
to lie on a cone whereby the order parameter 
space is a circle 5 , corresponding to the familiar 
[7(1) breaking. 

• Bose-Einstein condensates, which allow for a great 
variety of symmetry breaking if species of atoms 
are mixed. For example, consider two-species BEC 
(two different ultracold atomic gases in a trap) with 
order parameter fields <p a (x), (a = 1,2). The trap- 
independent part of the potential can be written as 

m 

V = Axd^l 2 - r,lf + A 2 (|^| 2 - r&) 2 +/3|^| 2 |^ 2 | 2 , (7) 

with tunable parameters, which is no more than 
J2J| rewritten for 50(4). For a single species (I s * 
term) we have an 0(2) symmetry that is totally 
broken, as in 4 He, and for two species an 0(2) x 
0(2) symmetry, broken to 0(2), with similarity to 
the breaking of the residual symmetry in high-T c 
superconductors . 



V. SYMMETRY BREAKING IN THE EARLY 
UNIVERSE 

In general, the overall patterns of symmetry breaking 
in the field theories that we believe describe the early 
universe are more complicated than those in condensed 
matter. In particular, the transitions about which we are 
most knowledgable (later in time) tend to be crossovers, 
and the ones we know least are likely to be complicated 
by virtue of their belonging to larger symmetry groups. 
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A. The Quark-Gluon System of Quantum 
Chromodynamics (QCD): 

Most known particles are strongly interacting 
(hadrons), of which the proton, neutron and pion are 
the lightest members. They are built from quarks 
(fermions) and gluons (vector bosons). Quarks are 
described by fields g/ iC , carrying two labels. The label 
c = 1,2,3 is that of a local SU(3) C symmetry, termed 
colour. In ordinary hadrons the gauge field gluons A a ^ c 
(a,c = 1,2,3), demanded by the local symmetry, form 
a massless octet (8 = 3 2 — l).The fact that there are 
are no free quarks in hadronic matter at zero chemical 
potential and temperature (particle accelerators) is 
explained if all known hadrons are colour singlets. The 
other label / = 1,2, ..,7V describes a global SU(N f ) 
flavour symmetry that is intrinsically broken. At Nj = 2 
this breaking is very small (electromagnetic), at Nf = 3 
it is still fairly small. For 3 < Nf < 6 the breaking 
is large. The resulting theory of quarks and gluons is 
known as Quantum Chromodynamics. 

The transition from a quark-gluon plasma to hadrons 
at low chemical potential (the early universe and heavy- 
ion collisions) is complicated by having two different as- 
pects. The first is the approximate breaking of chiral 
symmetry that is a consequence of the lightness of the 
common quarks, for which the gluons play a subsidiary 
role. The second is the confining/dcconfining nature 
of the gauge sector, for which the quarks are relatively 
unimportant. Although QCD can be tackled directly in 
some circumstances (through lattice gauge theory), it is 
informative to think of it as a modification of one of these 
extremes. More details can be found in the review article 
by Rajagopal 

• Low energy Chiral theory: Common hadrons 
are built from two quarks (termed up and down) 
which are approximately massless. One important 
idealisation is to take them massless in the first 
instance, md — m u — 0, to give a theory that is in- 
variant under independent transformations of the 
left and right hand components of the quark dou- 
blet ip = ipL+^R- Writing these as ipL — ► LipL and 
tpR — ► RipR, where L,R £ SU(2)f we would have 
degenerate parity doublets if the SU{2) l x SU {2) r 
chiral symmetry remains unbroken. This is not 
the case experimentally. Breaking the symmetry 
as SU(2) L x SU{2) R -> SU(2) L+R gives pions as 
the necessary Goldstone bosons. The pions are not 
massless (there are no known Goldstone bosons in 
particle physics) but they are anomalously light, 
requiring explicit (but small) symmetry breaking. 
This makes the transition a crossover. In fact, 
there is an additional U(1)a that is broken by in- 
stantons, without Goldstone bosons, that we shall 
not consider here. Extensions to Nf = 3 (requir- 
ing 'strange' quarks with m s — 0) with a more 
badly broken SU(3) L x SU(3) R -> SU(3) L+R will 



be made later. 

• Effective 0(4) sigma models: Still in the context 
of chiral theory, the chiral order parameter, a 2 x 2 
matrix M in flavour space, can be written in terms 
of a pion triplet 7?, with 7r + — tt\ + i7r2, etc, and 
a singlet a, as (^l a ^ R ) = M ab = aS ab + 7r.(f) ab . 
Low-energy pion physics is well represented by an 
effective 0(4) LaM , with potential V = \(tt 2 + 
a 2 - f 2 ) 2 (or an 0(4) NLaM with the constraint 
7? 2 + (T 2 = f 2 ) and explicit term linear in a to en- 
force TOtt ^ mass-breaking. In the absence of this, 
with manifold M. = S 4 of ground states the pions 
are the Goldstone modes (and the Higgs is the a). 
Note that 0(4) -> 0(3) w 0(3) x 0(3) -> 0(3) « 
517(2) x 017(2) -» SU{2). 

• Deconfinement If we take quark masses to be in- 
finitely large, we recover a pure gauge theory. In 
this case we can construct a non-local order pa- 
rameter (the Polyakov loop), whose sy mmetry- 
breaking characterises deconfinement [IJ- This is 
complicated by finite quark masses and is too sub- 
tle for us here. 

We should stress that there are many other ways to tackle 
the transitions, but they are even less relevant to our 
subsequent discussion. 

• Colour Superconductivity: At high baryon 
density (chemical potential) there is the possibility 
of Colour Superconductivity/ Superfluidity. This 
cannot be created in the laboratory but possibly ex- 
ists in neutron stars, which could have quark mat- 
ter cores. In cold dense quark matter the relevant 
degrees of freedom are those of quarks with mo- 
menta near to the Fermi surface. Quark attractions 
will lead to the creation of Cooper pairs. If we be- 
gin by assuming unbroken SU(3)f, with massless 
quarks m u — = m s = 0, Cooper pairs can- 
not be flavour singlets, and both colour and flavour 
symmetry is broken. The unbroken symmetry is 
(where L and R are chiral flavour) 

G = SU{Z) L x SU(3) R x SU(3) C x C7(1) B) 

where U(1)b counts baryons. The possibility ex- 
ists of colour-flavour locking (CFL), by condensates 

(4 , L a " l t'L b )3 x Ae Q/37 e° &7 , where a, a denote colour 
and flavour respectively. This assumes that the 3 
channel in 3x3 = 6 + 3 is the attractive channel. 

SU(3) L x SU(3) R x SU(3) C x U(1) B -> 

SU(3) L+R+c = SU(3) f+c . (8) 

Order parameters form a 3 x 3 complex matrix 
built from quark condensates (ipLipLjs and 
{4>riPr)3, as above, expressible as the nine Gold- 
stone bosons arising from the Chiral symmetry 
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breaking. There is an effective LoM realisation of 
this that we shall not pursue. 

Explicit symmetry breaking from mass terms is 
necessary. There are several possible phases as 
we further break the symmetry by breaking the 
SU(3)f quark mass degeneracy. On introducing a 
relatively massive strange quark, m s ^ m u = rrid, 
the K° is the pseudo-Goldstone boson, and readily 
forms a condensate. We have what is known as CFL 
+K°, in which SU(3) c+f SU(2) r X U(1) Y > -> 
U(l). 

The simplest scenario 01 has a doublet <f> of 

{ K + \ 

K°, K + charged bosons $ = ( ^ j f° r which 

V(«) = A(|$| 2 - 77 2 ) 2 (9) 
and M = 5 3 . 

If we include EM effects so that m s ^ m u ^ m c i 
we have the further breaking SU(3) c +f — > U(l) x 
U(l) — -> J7(l). The effect of the explicit mass break- 
ing terms is to recreate the potential V of Q for 
<h) = (K+, K°). 

B. Electroweak Transition/ Standard model: 

For systems of low baryon density, we have the un- 
ambiguous pattern of symmetry breaking of the Salam- 
Weinberg model, 

SU(3) C x SU(2)j x U(1) Y - SU(3) C x U(1) Q , 

where Q = I 3 + ^Y is the electromagnetic charge. This is 
replicated for each of the three families of quarks and lep- 
tons (corresponding to Nf = (1, 2), (3, 4), (5, 6)) respec- 
tively. Considering only the leptonic sector, the simplest 
scenario has a doublet $ of charged bosons 

$ = ^ for which V($) = A(|$| 2 - r] 2 ) 2 

of 0, after field relabelling. However, in this case V($>) 
is exact, with no further explicit symmetry breaking. 

On symmetry breaking we have the familiar three mas- 
sive vectors (W ± ,Z°) and one massless gauge field (the 
photon). The spin- 1 nature of the massive vectors has 
been bought at the expense of three of the field degrees of 
freedom of the $ doublet, to leave one massive real Higgs 
boson. For the assumed Higgs mass (m^ > 80Gey) the 
phase diagram shows a crossover. 

C. Grand Unified Theories: 

Ignoring supersymmetry, we assume the existence of 
symmetries G such that 

G -> Gsm = SU(3) c xSU(2) I xU(l) Y -» SU(3) c xU(l) Q . 



e.g E 6 -> 50(10) -► SU(5) -> SU(3) C x SU(2)i x U(l)y. 
In practice, we need an additional discrete symmetry to 
prevent too rapid a decay of the proton into leptons that 
is otherwise permitted by standard GUT symmetries. 
Groups like 50(10) contain a U(1)b-l gauge symmetry, 
which breaks to give an extended symmetry Gsm x Z2, 
forbidding fast proton decay. Motivated by the work 
of Witten 14] and others on superconducting strings , 
groups G containing U(l) x £7(1), with potentials (@J) and 
have been found useful, for reasons that we shall see 
later. 

However, once we attempt to accommodate supersym- 
metry at both a minimal level and beyond, the number 
of possibilities proliferates . 

VI. SIGNATURES OF TRANSITIONS 

The similarity between transitions in condensed matter 
and the early universe is interesting, but for it to be com- 
pelling we need a means of monitoring how transitions 
take place in condensed matter, so that we can check if 
similar assumptions about the early universe lead to our 
understanding it better. The way we do this is by look- 
ing for signatures of transitions from which we can infer 
the nature of the symmetry breaking. 

Consider the simplest case of a double-well poten- 
tial for a real scalar field <j>(x) We rewrite the potential 
V{4>) = |W> 2 - rff as 

V(6) — -(b 2 m 2 + —6 A + constant. 
2 4 

m = — A77 2 is negative. All this is at zero tempera- 
ture. As temperature increases and the field becomes a 
plasma m 2 becomes temperature dependent, increasing 
towards zero as m 2 (T) = to 2 (1 — T 2 /T 2 ) (in the simplest 
Ginzburg-Landau, or one- loop, approximation). The or- 
der parameter <j) — ((<j){x))), where the double average is 
both thermodynamic and quantum mechanical, satisfies 
4> 2 = rf(T) = ~m 2 (T)/X. At the critical temperature T c 
both (j) and m 2 (T) vanish. In relativistic quantum field 
theory T c ~ n, the symmetry breaking scale. For T > T c , 
m 2 (T) > and the Z 2 symmetry is restored. 

Let us now reverse the order, as in the early universe, 
and cool the field through T c from above, breaking the 
Z2 symmetry. How is it that symmetry is broken, given 
the 4> — <f> symmetry of the potential? The standard 
way, in field theory textbooks, is to bias the potential by 
introducing a spatially uniform external source j coupled 
to the field as j<p. For j ^ there is now a unique ground 
state to which the system relaxes. On taking j — > ±0 the 
system stays in this groundstate, with (f> — ±n according 
as the source is removed. 

In the early universe there is no such uniform bias. 
Instead, as we pass through T c there are local thermal 
fluctuations of the field that will drive it from the un- 
stable maximum at <f> — towards one or other of the 
minima. That is, for some x, (j>(x) is driven to +77, for 
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some it is driven to — T], and domains form in which <j> is 
correlated. In fact, since the transition is implemented in 
a finite time, causality (the finite speed at which the field 
can order itself) requires this domain structure. Once the 
system is behaving classically, the boundaries of these do- 
mains are 'walls', in which the field passes from <f> = —r\ 
to (j> = r\ according to the classical equation 5S/5<p = 0. 
Such a 'domain wall' in the x — y plane, with profile 
4>{z) = 77tanh(|m|z/\/2), has thickness £o ~ l^l" 1 and 
energy per unit area (surface tension) a ~ |m|?7 2 . We 
will know that the transition has taken place by the exis- 
tence of these walls. Already we have the powerful result 
that, if we do not adopt an inflationary scenario to dilute 
wall formation, then domain walls produced at the GUT 
scale are so energetic that they would close the universe, 
in contradiction to the evidence. Although a closed do- 
main wall can contract and disappear, it is not possible 
to remove part of a wall, and a single wall requires an 
infinite field reordering to remove. Domain walls are the 
simplest entities to show this topological stability, and 
they provide the simplest example of a topological defect. 

More generally, as we shall see, topological defects are 
twists, or knots, or other deformation of the fields that 
carry a topological charge that prevents them unwinding, 
or untwisting (but not annihilating). Whether defects 
exist depends on the detailed nature of the symmetry 
breaking. If they do exist, individual defects produced in 
the early universe can survive into the present era as 'fos- 
sils'. Observing them would provide concrete evidence 
for transitions having taken place and, if they could be 
studied, provide information on the symmetry breaking. 
We can check this out directly with condensed matter 
systems that permit defects. Further, we can make pre- 
dictions for defect formation that can be checked in con- 
densed matter, and which have implications for the early 
universe. 

Of course, for different patterns of symmetry breaking 
the frustration in the order parameter fields as they link 
domains formed after the transition may not be resolved 
with topological defects. In such cases the twists in the 
fields can unwind at no great energy cost and, at times 
well after the transition, there is no direct signal for it 
having taken place. In practice, models for the early 
universe almost inevitably produce defects of one form 
or another, and the problem is more one of not seeing 
them. 



VII. TOPOLOGICAL DEFECTS 

A. Vortices and Strings 

The most common topological defect in condensed 
matter physics and the early universe is the string, or 
vortex, seen most simply in the breaking of a U(l) sym- 
metry. 

Consider a global theory of a complex scalar field 
4>(x) — p(x)e l9{ - x \ with effective potential V(<p,T) — 



j(\(f>\ 2 — rj 2 (T)) 2 , where rf(T), vanishing at the critical 
temperature T c , is given above. As we cool through T c 
the field, originally concentrated around the (then) sta- 
ble minimum at <fi = 0, begins to explore the possible 
groundstates with p = rj. With no bias to make <j> fluc- 
tuate away from the unstable maximum at <fr = with 
any particular phase, 9{x) will vary from point to point, 
subject to continuity. 

Let us take a closed path C in space. Once the tran- 
sition has been completed, \<f>{x)\ w r\ for x e C. If we 
make a complete circuit in C then <fi(x) executes a closed 
path in S 1 , the circle of ground states. Continuity of 
4>(x) requires that the change A8 in 9{x) along the path 
is an integer multiple of 2-k, A8 = 2nn. We can find 
classical solutions of infinite extension along the z-axis 
of the form <p n (x) = p(r)e m9 , where r 2 = x 2 + y 2 , and 
6 = tan _1 (a;/y) is the azimuthal angle, that demonstrate 
this. In 4 He they are vortices with quantised vorticity, 
of winding number n. It happens that, if n > 1, it is 
energetically advantageous for the vortex to break into n 
vortices of winding number unity. These are the stable 
topological defects of the model, being unbreakable. If 
we shrink the loop C to the origin r = 0, we must have 
Pn(0) = for 4>(x) to be defined. That is, the core of 
the vortex is the false groundstate, just as the core of do- 
main wall was. In low-T c superconductors the coupling of 
the electromagnetic and Cooper pair fields makes them 
Abrikosov vortices with quantised magnetic flex. As hap- 
pens for 4 He, the Cooper pair field (j> — at the vortex 
core. The thickness of the false vacuum (cf> w 0) core 
is the London length, the thickness of the flux tube the 
Meissner length. In the early universe, with no Gold- 
stone bosons, in general vortices will be of the Abrikosov 
type. 

Suppose we have collection of vortices of the type 
above. Consider two closed paths C\ and C2 with a point 
in common, in which the field takes winding numbers 
m and ri2- We can trivially combine them into a closed 
path C\ + C2 with winding number m + n 2 - That is, the 
mappings from closed paths Cj, which we can take to be 
circles S 1 into the set of ground states M = S 1 form the 
additive group of integers Z. The group Z is termed the 
Fundamental Group or the First Homotopy Group of S 1 , 
written H^S 1 ). 

More generally, a transition will produce topologically 
stable vortices if Hi(M) = Ii\{G/H) is non-trivial. For 
the condensed matter systems that we have discussed 
above, we also find topological Z vortices in single-species 
BEC and liquid crystals at the interface of an isotropic- 
nematic transition. For our broken 50(5) high-T c super- 
conductor the situation is more complex, although vor- 
tices have winding number n € Z . As we get close to 
critical doping, and SO (5) symmetry is restored (prior 
to coupling to elcctromagnctism) the stable vortices have 
non-trivial antifcrromagnctic cores in which the order pa- 
rameter 77 7^ 0, although necessarily vanishes there. 
As we achieve critical doping the antiferromagnetic core 
expands to destroy the vortex, since IT (.J? 4 ) is trivial, not 
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permitting topological charge. 

For superfluid 3 He the situation is even more com- 
plicated. In the ,4-phase, ^(S 2 x SO(3)/Z 2 ) = Z 4 , 
addition modulo 4. There are four types of vortices, 
with winding numbers n = 0, ±1/2, ±1 (mod 2). De- 
pending on the parameter values, it is again possible 
to have vortices with non-trivial cores. In the i?-phase, 
III (S 1 x SO (3) ) = ZxZ%, and there are two types of vor- 
tices with arbitrary winding number. Not all vortices are 
normal (with all components of A a i vanishing at their 
cores). It is possible for S-phase vortices to have A-phase 
cores. 

For cosmology, topological strings (vortices) appear in 
colour superconductivity and superfluidity. In supersym- 
metric GUT transitions which solve the monopole prob- 
lem, lead to baryogenesis after inflation and constrain 
proton decay, the vast majority of models lead to topo- 
logical (cosmic) strings p~5j . 

We should not ignore so-called embedded strings, in 
which, when G — > H, there are subgroups G", W of G, H 
respectively, for which TL\(G' /H') is not trivial. Such 
strings occur in the 0(4) sigma model of low-energy pi- 
ons and in electroweak breaking 5/7(2) x U(l) — ► U{1). 
However, such strings have no topological winding num- 
ber and whether they are energetically classically unsta- 
ble or not depends upon the values of the parameters 
(masses, coupling constants). Even if they are the role of 
quantum fluctuations is not clear. In the important case 
of electroweak breaking the strings are unstable to begin 
with, and can only have transitory effects. 

Loops of simple relativistic strings contract and turn 
their energy into particle production. Although vortex 
loops in A He can be stabilised by Magnus forces, in gen- 
eral loops conract. However, many theories (most simply 
of the U(l) x U(l) — > U(l) type) can produce strings or 
vortices that have non-trivial cores with angular momen- 
tum, which can stabilize loops (vortons) Such a 
mechanism is not possible in 3 He on energetic grounds, 
although it is possible in 2-component BEC, where it is 
possible that an example has been seen. The main case 
of theoretical interest is that of CFL+if quark matter, 
because of its relevance to neutron stars. 



B. Monopoles and more 

Consider an 50(3) sigma model with 3-vector field 
$(x) and potential V{$) = A/4(0 2 - r] 2 (T)) 2 as it cools 
through its critical temperature at which rj{T c ) = in 
a finite time. Causality requires the direction of <p to 
be uncorrelated at large distances. This frustration is 
resolved by classical 'knots' in the field, of the form 
4>{x) = p(r)x/r, where r is the radial distance from the 
centre of the knot. These are the topological monopoles 
of the model and cannot be 'unwound'. If we enclose one 
by a sphere S 2 , we see that they are non-trivial map- 
pings of S 2 onto the manifold of ground states M. = S 2 . 
These mappings again form a group, the Second Ho- 



motopy group n 2 (A4). In fact, n 2 (5 2 ) = Z, and the 
monopole solutions are characterised by integer winding 
number, for which our example (|n| = 1) is stable. 

More generally, we have monopoles whenever n 2 (A4) 
is non-trivial. In condensed matter physics we find 
monopoles in the A-phase of superfluid 3 He, for which 
n 2 (5 2 x 50(3)/Z 2 ) = z. 

Monopoles arise inevitably in GUT transitions. This 
is a consequence of the fact that, if G is connected and 
semisimple (i.e. IIi(G) = 1), thenn„(G/iJ) = n„_i(If). 
Thus, if Hi(H) 7^ 1, we have monopoles. Since H neces- 
sarily contains U(1)q, Hi(H) ^ 1. Just as domain walls 
produced in a non-inflationary universe would close it, so 
would monopoles. It was to dilute the monopole density 
to an acceptable level that Guth originally introduced 
inflation. This does not mean that all defects are di- 
luted. There is no difficulty in developing hybrid models 
in which monopoles are diluted but strings survive |15|. 

All the above is only the simplest possibilities. With 
suitable symmetry breaking we can have strings ending 
in monopoles, domain walls bounded by strings, and so 
on. 



C. The observation of defects 

Let us postpone estimating how many defects we might 
see at a transition until later, and ask the simpler ques- 
tion as to whether we see examples of these defects in 
real life. All the simple defects permissable in condensed 
matter systems are seen in the laboratory, as are many 
more complicated combinations of defects. See articles 

In cosmology the situation is the converse, with no un- 
ambiguous sightings of monopoles or strings (vortices). 
Cosmic strings can be best observed by gravitational 
lensing, in which no central object is seen. The met- 
ric around a cosmic string is conical, with a defect angle 
5 = 8tt(T c /M P i) 2 ~ 10~ 5 . Recently a plausible galaxy 
string lensing candidate was observed |17|. but in the 
absence of further lensing nearby, the existence of a cos- 
mic string is in doubt. Cosmic strings are, in principle, 
a potential source of Extra High Energy Cosmic Rays 
(EHECR) ^^|. However, simple mechanisms like loop 
decay, cusp radiation and string interconnection give a 
flux many orders of magnitude (> 10) too small. There 
are many other variants, among which vorton loops and 
loop s of superconducting string are the more promising 
|l9| . In the latter case, the superconducting current in- 
creases until the loop disintegrates, producing EHECR. 

The most promising source of quark-hadron transitions 
is the interior of neutron stars |lll |20j ■ The quark-gluon 
condensate is highly compressible, compared to nuclear 
matter. Arguably this is seen in millisecond pulsars. As 
they spin-up via accretion the core-density decreases and 
the core turns to hadrons. As the quark-matter spins 
it inhibits an increase in the pulsar frequency u) because 
of an increase in the moment of inertia. Detailed calcu- 
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lation suggests a frequency (uj ~ 300Hz) at which the 
pulsar distribution should peak, and this is observed ex- 
perimentally |2lj . 

The vortices and vortons that are produced at the tran- 
sition provide candidates for the observed glitches in pul- 
sar frequencies (Slu/uj < 10~ 6 ). The suggestion is that 
these are understood as a consequence of the interac- 
tions between the rigid crust of the neutron star and vor- 
tices/vortons in the neutron superfluid interior as they 
try to move outwards. 



where a = jis. This is very large on the scale of cold 
defects. We term a the Zurck-Kibblc (ZK) characteristic 
index. In the mean-field approximation a = 1/4 typically 
for condensed matter, and a = 1/3 for QFT. 

We see that the same arguments apply to crossovers 
p^j , provided they are weak enough that the correlation 
length is larger than the causal horizon at time t, even 
though it does not diverge. We assume this to be the 
case. 



VIII. THE DYNAMICS OF DEFECT 
PRODUCTION 

A. Causal Bounds 

Because phase transitions take place in a finite time, 
causality guarantees that correlation lengths remain fi- 
nite. If the symmetry breaking permits, defects arise so 
as to mediate the correlated regions with different ground 
states. They provide an excellent experimental signature 
for the way in which transitions are implemented. 

For condensed matter systems, Zurek|22l 23] suggested 
that causality alone is sufficient to determine the initial 
density of defects arising in a continuous transition. In 
this he paralleled proposals made by Kibble[24| in the 
context of quantum field theory models of the early uni- 
verse. 

Consider a system with critical temperature T c , cooled 
through that temperature so that, if T(t) is the temper- 
ature at time t, then T(0) = T c . T(0) = T c /tq defines 
the quench time tq. One formulation (see |23jp of the 
Zurek/Kibble causal bound is the following: Suppose, 
at time t, that c(t) = c(T(t)) is the maximum speed at 
which the order parameter can respond to a change in 
the environment. Then, at time t the causal horizon has 
diameter 



h{t) 



ds c(s). 



(10) 



Let £ad(t) = £ad(T(i)) be the adiabatic healing length 
(correlation length), diverging at time t = (T = T c ). 
The first time that defects can appear is at time t, when 
the causal horizon becomes big enough to enclose one of 
them, £ od (i) = h(t). 

As a direct consequence, t has the form i — Tq 7 Tq . 
The scaling exponent 7 depends upon the system. For 
dissipative condensed matter, with critical slowing down 
(c(0) = 0), r <C tq is the relaxation time of the longest 
wavelength modes, for QFT (with c a constant) the time 
for light to traverse a Compton wavelength. As a result, 
tq > t > tq. If £ ad (t) ~ £o(*/tq)-" for t ~ 0, where 
£0 = 0{^ a d{T — 0)), then the initial domain size and 
defect separation is predicted to be 



B. Unstable modes 

In the several years since these simple bounds were first 
proposed we have acquired a much better understanding 
of the way in which transitions occur. These does not 
mean that these bounds have lost their relevance, but 
that they need to be qualified. 

Any dynamical equations for the onset of a continuous 
transition will embody causality, by definition. However, 
the transition cannot be said to have happened before 
the order parameter has achieved its equilibrium value 
\(j)\ 2 = rj 2 , (in our notation for potentials). If (...)* denotes 
ensemble averaging at time t then a lower bound on the 
first time from which we can start counting defects is 
t = **, for which {\(j)\ 2 )t = rj 2 . 

For many condensed matter systems and for quantum 
fields, the way in which (\4>\ 2 )t builds up to its final value 
is by the growth of the amplitudes of the unstable long- 
wavelength modes, which are unstable because of the up- 
turned parabolic free energy density at initial times. The 
time t* is, crudely, the time for these modes to roll from 
the top of the hill to the groundstates at the bottom, 
related to the spinodal time. 

The equations that control this ordering through in- 
stabilities differ for different systems. 



C. Dissipative systems: 

In the vicinity of T c , dissipative systems are often 
well approximated by the time-dependent Ginzburg- 
Landau equations (TDGL). For a global theory these are 
Langevin equations of the form 



r dt 



SF 



(12) 



T 



(11) 



where rj a is Gaussian thermal noise, satisfying 
(»7«(x,t)»fc(y,f )) = 2S ab T{t)T5(x - y)6{t - f). This is 
a crude approximation for A He, and a simplified form of 
a realistic description of 3 He pj, and an important part 
of the description of low-T c superconductors. 

This seems a very different picture from that originally 
proposed by Zurek. A priori, t* is not related directly 
to the causal i, but it is not difficult to see why they 
might be comparable. Unstable modes grow exponen- 
tially fast. As long as dimensional analysis makes i the 
natural unit in which to measure time, any exponentially 
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growing term will achieve values that are not exponen- 
tially large at times t — 0(t). This can be checked explic- 
itly for fast quenches on retaining only the quadratic part 



of F, assuming that growth in 



is then checked by 



the (linearised) back-reaction. This is justified by seeing 
that the growth in long wavelength modes is largely com- 
pleted in the linear regime, when the field components (f> a 
are effectively independent [26|, |27j • 



D. Relativistic QFT: 

For QFT the situation is rather different. In the previ- 
ous section, instead of working with the TDLG equation, 
we could have worked with the equivalent Fokker-Planck 
equation for the probability pt[$] that, at time t > 0, 
the measurement of <fi will give the function ^(x). When 
solving the dynamical equations for a hot quantum field 
it is convenient to work with probabilities from the start. 

The probability pt [$] that, at time t, the measurement 
of <fi will give the value $ is p t [$] = | ^> | 2 , where is the 
state-functional with the specified initial condition (e.g. 
Boltzmann distributed). In the language of path inte- 
grals, pt[<&] can be written as an integral for fields with 
increasing time followed by an integral for fields with 
decreasing time (\f r *). The time contours can be joined to 
give a closed time-path and non-equilibrium calculations 
in QFT are termed closed time-path calculations. 

In practice, there is no need to calculate Pt[&\ directly. 
If (|$(a;)| 2 )t measures the growth of field modes as an 
ensemble average with respect to pt [ < I > ] , then 

(\<f>(x)\ 2 ) t = {\4,(x,t)\ 2 ) 

for Wightman fields cj)(x,t), subject to the thermal 
boundary conditions at the initial time before the transi- 
tion is implemented. Mode analysis is all but impossible 
outside self-consistent linearisation of the back-reaction 
pjjj . Within the self-consistent linear regime (which is 
the best that can be done numerically) the mode equa- 
tions that determine (|</>| 2 )t are now the classical second- 
order equations 



d 2 



SF 



dt 2 



(13) 



Yet again, this looks a rather different picture from 
that originally proposed by Kibble, in which the system 
froze in before the transition was effected. Nonetheless, 
although t* is not related directly to the causal t, they 
are again comparable, essentially because unstable modes 
grow exponentially fast, and t has the correct engineering 
dimensions. 

There is one major difference between condensed mat- 
ter systems and the early universe that is not addressed 
in equations like (|13fl . This is that, whereas the effec- 
tive Ginzburg-Landau field theories for condensed matter 
systems are complete, the early universe contains many 
(probably most) fields about which we know nothing, and 



which are ignored in JT3J. Only for QCD to we have a 
full tally of the relevant degrees of freedom of the sys- 
tem. The effect of this environment is to make the quan- 
tum field theory behave classically (show decoherence) 
[2^. l30| . The easiest way to show this is by seeing how 
the density matrix (whose diagonal elements are pt[$]) 
becomes diagonal as a consequence of the growth of the 
unstable modes. In parallel, the master equation for the 
Wigner functional plays the role of a Fokker-Planck equa- 
tion, whose Langevin counterpart for the semiclassical 
field is a variant of (|12|l . including multiplicative noise. 
However, none of this interferes strongly with defining t* 
from (Q3J). 

What is important is that, by the time that defects are 
produced they behave like classical entities |scj . 



E. Other systems 



Not all systems behave as the above. Bose-Einstein 
condensates satisfy the Gross-Pitaevski equation ^(| 



SF 



dt 



(14) 



Although causality still leads to the Zurek bounds, 
transitions do not evolve because of the exponential 
growth of unstable modes. A similar situation exists for 
those transitions that are driven by changes in chemical 
potentials (as in high-T c superconductors). 



F. Additional mechanisms for gauge theories 

All the mechanisms above are well-suited for global 
symmetry-breaking, but for local symmetry-breaking in 
the presence of gauge fields there is an additional mech- 
anism for the production of defects. This mechanism, 
observed by Hindmarsh and Rajantie (3 1| is discussed in 
detail in Arttu Rajantie's contribution to these proceed- 
ings, and will only be considered briefly here. 

For the most relevant case of superconductors (both 
high and low temperature) the idea is simple, even if 
the execution is difficult. Above the critical tempera- 
ture there are thermal fluctuations in the electromagnetic 
field. As the system cools through T c the short wave- 
length modes of the field stay in equilibrium, whereas 
long wavelength modes drop out of equilibrium and freeze 
in. These frozen modes are the source of correlated flux 
that will be measured along with the flux of conventional 
Abrikosov vortices produced from the causal arguments 
above. The effect is not simply additive but will, in gen- 
eral, lead to a greater variance in the spontaneously pro- 
duced flux passing through a surface. 
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G. Defects as zeroes 

So far we have only been discussing the time t* it takes 
for the transition to complete itself, and comparing this 
to the causal time t. The important step is to derive the 
defect separation £* at the time of defect production and 
compare it to £ of l(TT|) . 

To do this, we identify simple defects, with normal 
(false vacuum) cores, by the zeroes of the order parameter 
fields. At early times, when field fluctuations are strong, 
field zeroes do not have the energy profiles to be identified 
with defects, whose masses, tensions, etc, are 0(rf) non- 
perturbatively large. However, by t*, the energy profiles 
are qualitatively correct, and counting zeroes provides a 
reliable estimate of defect numbers. [In the case of QFT, 
this requires that defects can be considered classically] 

In the linear regime the separation of zeroes is given 
|32j | simply in turns of the derivatives of the two-point 
correlation function (4> a (x, t)4>b{0, t)). In this approxi- 
mation we find that, just as t* ~ t (up to logarithmic 
corrections), then so is £* « 

We have observed that many systems have defects with 
complicated cores. Simple numerical simulations show 
that, provided the non-normal components of the cores 
are not too la rge, the scaling laws of the causal analysis 
are preserved [33| . 

There is a further important comment on thermal 
fluctuations. Initially |34| . it was thought that domain 
size after a transition was determined by the correlation 
length at the Ginzburg temperature T c . In the causal 
arguments above the Ginzburg regime plays no role, but 
this reappears in the analysis of unstable modes, where 
thermal fluctuations on small scales can make the pro- 
duction of defects difficult [35| if the quench is too slow. 



IX. EXPERIMENTAL CONFIRMATION 

To date, several condensed matter experiments had 
been performed to test ifTTl) : 

• Superfluid 3 He - H^. The two 

experiments [3^, |37| on 3 He — B rely on the fact 
that, when superfluid 3 He — B is bombarded with 
slow neutrons, n + 3 He — > p + 3 H + 760keV. The 
energy released in such a collision leads to a hot 
spot in the superfluid, with temperature T > T c , 
which when cooled by its environment, leaves be- 
hind a tangle of vortices (the topological defects in 
this system), tq is fixed by the nuclear process that 
breaks up the 3 He atom. With only a single data 
point conflating both normalisation and a it is not 
possible to confirm the predicted value a — 1/4. 
However, both experiments are highly compatible 
with ||TTJ|. 

• Superfluid In principle, the two 4 if e 
experiments |38l l39| . which use a pressure quench, 



allow for a more complete test. Yet again, vortices 
are the relevant defects. In practice, the most re- 
liable experiment j3^ sees no vortices. This is not 
necessarily a sign of contradiction in that it has 
been suggested 35] that the vortices decay too fast 
to be seen in this case. This is irrespective^^ °f 
whether high thermal fluctuations within the wide 
Ginzburg regime of 4 He would lead to somewhat 
different predictions. In this context, the vortices 
seen in an earlier 4 He experiment |38j were most 
likely an artefact of the experimental setup. 

• High temperature superconductors 
(HTSC)Hl Eg. The two experiments|i3, 
on HTSC measure the total flux of Abrikosov 
vortices through a surface. The vortex separation 
of (fTT|> can be converted into a prediction for the 
variance. In the first experiment no flux was seen, 
despite the phase separation that leads to the 
result being demonstrated elsewhere |43j. However, 
on increasing the quench rate by several orders of 
magnitude spontaneous flux was produced [42^ . 
but with low efficiency. The observed flux is 
compatible with the Zurek prediction of (| 1 1|) . once 
additional assumptions for converting total flux 
into net flux are taken into account. Although 
the prediction l|ll|) has not taken gauge fields 
into account, the effect of electromagnetic field 
fluctuations is expected to be small [44^. 

• Annular Josephson Tunnel Junctions. The 

defects of a lineaer JTJ are fluxons, the sine- 
Gordon model kinks in the field <f> = <f>\ — 4>2, 
the difference in the phases of the complex or- 
der parameter fields in the separate superconduc- 
tors. Again two experiments have been performed 
(one currently in progress, with data as yet un- 
analysed). The first experiment was one in 
which, by counting fluxons on varying quench time 
tq, it was possible to compare a with its theo- 
retical value (as well as confirming overall scale). 
The secondary mechanism of [3lJ is n °t relevant 
here. Agreement with both exponent and mag- 
nitude are good, although there is scatter, with 
a = 0.27 ± 0.05, in comparison to the theoretical 
value of a = 0.25. The new experiment currently 
being performed will permit quench times four or- 
ders of magnitude faster than before, and thereby 
an order of magnitude more fluxons. 

• Other direct experiments. Two subsequent ex- 
periments have permitted varying quench rates and 
so an estimate for a. The first |46j involves the 
Benard-Marangoni conduction-convection transi- 
tion, in which a homogeneous conduction state is 
broken into an hexagonal array of convection lines 
on heating. The defects here are not associated 
with the line zeroes of an order parameter field, 
and the viscosity-dependent a does not match the 
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ZK prediction, possibly for that reason. The second 
experiment |47j is carried out in a non-linear opti- 
cal system, with complex beam-phase the order pa- 
rameter. Increasing the light intensity (the control 
parameter in this case) leads to pattern formation 
(defects) at a critical value. The predicted a = 1/4 
is recovered to good accuracy as a exp = 0.25 ±0.02. 

• Related experiments. Related experiments in- 
clude measurements of defect production in con- 
tinuous transitions of planar liquid crystals 0, \^ , 
measurements of flux production in dilute Joseph- 
son Junction annuli 49] , and measurements at first 
order transitions in radioactive low-T c supercon- 
ductors ■ The first provides a check on the rela- 
tionships between net and gross topological charge, 
the second on the role of electromagnetic field fluc- 
tuations at a transition |44| , relevant for the high- 
T c experiments. At the moment the implications of 
the third are unclear, given the different nature of 
first order transitions. 



X. CONCLUSIONS 

We have seen that the patterns of symmetry break- 
ing in condensed matter systems and the early universe 



are sufficiently rich and sufficiently similar to pursue the 
comparison. 

In almost all cases the transitions are signalled by the 
production of defects, usually topological vortices, vor- 
tons or monopoles. The observation of defects is therefore 
a demonstration that the transitions have taken place. 
Defect production in the early universe is complicated 
by the assumed inflationary period, and there are only 
tantalising glimpses of possible candidates. The situation 
is very different for condensed matter systems, where de- 
fects are readily observed. 

In this latter case we can do more, and examine the 
details of defect production. A comparison of the ex- 
perimental results with the causal constraints on defect 
production proposed by Zurek shows the strength and 
limitations of the causal bounds and informs the similar 
arguments on the role of causality in defect production 
in the early universe, as originally mooted by Kibble. 
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